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Abstract 

We consider 2+1 flavor Wilson Chiral Perturbation Theory including the lattice spacing contributions of 0(a 2 ). We adopt 
a power counting appropriate for the unquenched lattice simulations carried out by the CP-PACS/JLQCD collaboration and 
compute the pseudo scalar meson masses to one loop. These expression are required to perform the chiral extrapolation of the 
CP-PACS/JLQCD lattice data. 
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I. INTRODUCTION 



The limitations of the quenched approximation in numerical lattice QCD simulations is by now well established. For 
example, the light hadron mass spectrum calculated by the CP-PACS collaboration pj deviates from the experimen- 
tally measured values by about 10 percent. Even though the quenching error is different for different observables, one 
must assume the quenching error to be of the same order for other quantities as well. Once the effect of dynamical up 
and down quarks is included, the quenching error is significantly reduced and the discrepancy between the numerically 
calculated and the experimentally measured values is much smaller compared with the quenched results Q. Still, 
ignoring the effect of a dynamical strange quark in these unquenched 2 flavor simulations leads to an uncertainty, 
which is expected to be non negligible. Only simulations with a dynamical strange quark will provide numerical 
results, which can be compared with experiment with confidence. 

In order to eliminate the remaining source of quenching error the CP-PACS and JLQCD collaborations have been 
carrying out unquenched 2+1 flavor simulations. A RG-improvcd gauge action and an 0(a) improved Wilson quark 
action have been adopted. The size of the lattice is modest (L ~ 2 fm) and simulations at three different lattice 
spacings (a ~ 0.7 fm, 1.0 fm, 1.22 fm) are planned so that the continuum limit can be taken. Five different masses 
for the degenerate up and down type quarks are simulated leading to pseudo scalar meson masses in the range 
mps/my ~ 0.62 — 0.78. The physical strange quark mass lies between the two simulated strange quark masses and 
can therefore be reached by an interpolation. More details and the status of these simulations have been recently 
summarized in Ref. Q- 

The masses for the up and down quarks are rather heavy and an extrapolation to their physical values is required. 
The functional forms for the extrapolation are usually given by Chiral Perturbation Theory (ChPT). In its standard 
form U,[5| the expressions derived in ChPT can be used after the continuum limit of the lattice data has been taken. 
However, for various reasons it is advantageous to perform the chiral extrapolation before the continuum limit. In this 
case ChPT needs to be formulated for lattice QCD at non-zero lattice spacing a. The main idea how this can be done 
was proposed in Ref. 0, Q ■ Since then many observables have been calculated at one loop order (for an overview see 
Ref. ||)- For lattice theories with Wilson fermions, however, all results were derived for unquenched 2 flavor Wilson 
ChPT (WChPT) and the chiral expressions are therefore not applicable for the CP-PACS/JLQCD simulations. 

This is the first paper in a series where we provide the one loop expressions of 2+1 flavor WChPT for a variety 
of observables, which will be measured by the CP-PACS/JLQCD collaboration. Here we present the results for the 
simplest observables, the pseudo scalar meson masses. The second paper is devoted to the vector meson masses and 
the third to the pseudo scalar decay constants and axial vector Ward identity quark mass 0, El . In our calculations 
we include the lattice spacing contributions through 0(a 2 ) and adopt various power countings. Even though we have 
primarily the CP-PACS/JLQCD simulations in mind, our expressions are equally useful for other unquenched 2+1 
flavor simulations employing Wilson fermions. 

There is no fundamental difficulty in applying the framework of WChPT to 2+1 flavors, the main difference to the 
2 flavor results is just the increased complexity of the final results. Since we follow the standard strategy of WChPT 
we will be brief in presenting our results. In section ^ we explain the power counting which we assume and present 
the chiral Lagrangian up to next-to-leading order. The calculation of the pseudo scalar masses from this Lagrangian 
is straightforward and we summarize our final results in section ITTT1 Many technical details and some intermediate 
results are collected in appendices lAl and iBl 

II. THE CHIRAL EFFECTIVE LAGRANGIAN 
A. The order counting 

In continuum chiral perturbation theory (ChPT), M — m or p 2 is the expansion parameter, where m is the quark 
mass and p is the momentum of the pseudo scalar meson. Since chiral symmetry is explicitly broken in lattice QCD 
with Wilson-type quarks, corrections due to the non-zero lattice spacing a are non-negligible. The construction of the 
chiral effective theory for Lattice QCD with Wilson fermions, so-called Wilson ChPT (WChPT), has become standard 
by now. From a conceptional point of view there is nothing new in applying the familiar techniques 0, ITU IT3 . Il3l | to 
2+1 flavor lattice QCD. The only non-trivial choice one has to make is the order counting one adopts in the chiral 
expansion, which we are going to explain in this section. 1 



1 For notational simplicity we assume N degenerate quark masses in the following discussion. 
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The leading order (LO) chiral Lagrangian of the WChPT is constructed from O(M) terms and the 0(a) term. 
Since the 0(a) term has the same chiral structure as the mass term, the LO Lagrangian of the WChPT assumes the 
same form as the one of continuum ChPT\provided one performs the replacement m — > rh = m + cia with c± being a 
combination of two low-energy constants |6llll|. Based on this LO Lagrangian, however, the pion becomes a tachyon 
for rh = m + cia < since the pion mass squared m T is given by m„ = 2B(m + cia) at tree level. This is in contrast to 
continuum ChPT where chiral symmetry dictates that the pion mass is given by m 2 = £?|m|. This already indicates 
that the first non-trivial modification to continuum ChPT starts at 0(a 2 ). Indeed, including the 0( a 2 ) term in the 
LO chiral Lagrangian removes the unphysical tachyon from the theory, as has been shown in Ref. [12( for the 2 flavor 
theory. We therefore conclude that, for the consistency of the theory, the 0(a 2 ) term should be included in the LO 
Lagrangian in the WChPT. 

Although the 0(a) term is superficially larger than the 0(a 2 ) term, the former term is irrelevant since it is absorbed 
in the definition of the shifted mass rh. The 0(a 2 ) correction, on the other hand, becomes important in the regime 
where rh /Aqcd — Aq CD a 2 , even though the Aq CD a 2 correction is much smaller than 1 in general. 

Suppose we consider the 0(M) = 0(m,p 2 ) and the 0(a 2 ) term as LO terms. Then the terms of 0(M 2 ,Ma 2 ,a 4 ) 
can be regarded as next to leading order (NLO), since the loop expansion in WChPT increases in units of M. We 
remark that the 0(a ) term is not as relevant as the 0(M 2 ,Ma 2 ) terms for our final results for the pseudo scalar 
masses. We will need the tree level contribution of the 0(a 4 ) term only to satisfy rrip S NLO = at m,p S LO = 0, 
where mps.NLO is the pseudo scalar meson mass at NLO. This means that there is essentially no unknown low energy 
constant associated with the 0(a 4 ) correction. 2 

The proper order counting of the O(Afa) term is more subtle than for the previously discussed terms. 3 Depending 
on the size of the O(Ma) contributions we may include it at leading order where it subsequently enters the chiral 
logs, or we treat it as a NLO term and include it at tree level only. The size of this term is indeed expected to vary 
significantly, depending on details of the action of the underlying lattice theory. The O(Ma) term contains one power 
of a and stems from the Pauli term in the Symanzik's effective action |l6l Il7l HS|. which is used in an intermediate 
step to match the lattice to the chiral effective theory. Conse quen tly, the O(Ma) contribution in the chiral Lagrangian 
is directly proportional to the coefficient of the Pauli term [ll|, and the size of this coefficient is much smaller for 
improved theories with a clover term in the action than for standard Wilson fermions. For fully non-pcrturbativcly 
improved theories the coefficient is equal to zero and no O(Ma) term appears in the chiral Lagrangian. 4 

Since we have no a priori knowledge about the size of the O(Ma) term we will be as general as possible and present 
our results for three different order countings of this term. We first keep the O(Ma) term at LO where it gives a 
contribution to the chiral log corrections. If the O(Ma) term is assumed to be much smaller than the other LO 
terms we can easily expand our results. Performing this expansion is equivalent to doing the calculation with keeping 
the O(Ma) term at NLO. Finally we can set this term to zero in order to obtain the results for non-perturbatively 
improved theories. 



2 This may seem odd within WChPT, but it is a simple consequence of the fact that the pseudo scalar meson becomes massless in lattice 
QCD at the critical quark mass. In other words, the 0(a 4 ) term merely results in an additional shift in the cri tical quark mass. An 
implicit assumption we make here is the presence of a phase where flavor and parity is spontaneously broken |14| . 

3 The 0(a 3 ) term is unproblcmatic since the arguments we gave for the 0(a 4 ) contribution also apply for the 0(a 3 ) term. 

4 Strictly speaking this only holds if the chiral Lagrangian is parameterized in terms of renormalized quark masses which absorb some of 
the 0(a) cut-off effects. Particular terms linear in a will appear if other choices are used (see section fill CI . 
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B. Leading order Lagrangian 



According to the discussion in the previous section, we use the following LO Lagrangian, which consists of terms 
of 0(M, a 2 , Ma): 

+ |! [c 2 (S + St) 2 + C ~ 2 ((E + Et ) 2 ) + C4(s _ E t)2] 
16 

f 2 

+ J — [c (E + S f - 2) (^s^st) + £ <(S + s f - 2)a„Ea„Et)] 

f 2 i? 

+ -IT [ 2c s( E + S + )(M g E + E+M,) + c 3 ((E + Et)(M,E + E+M ? )}] 

o 

+ ^ [ 2c 5(S - St)(M g E - Et M ,) + c 5 ((E - Et)(M g E - E+M,))] , (1) 



where (X) = tr X , / is the pscudo scalar meson decay constant, 



cxp 



(2) 



is an element of SU(3) with 7r a being the pseudo scalar meson fields. The SU(3) generators T a are normalized 
according to tr T a T h = ^8 a b- The first and the second term in the first line are the standard 0(p 2 ) and O(m) terms 
0, respectively. The second line comprises the 0(a 2 ) terms [HQ. The last three lines contain the 0(p 2 a) and 
O(ma) contributions [Til ]. 

For notational simplicity only we use a different notation for the low e nerg y constants associated with the non-zero 
lattice spacing (the c and c's) compared to the notation in Refs. [ill [l3l| . In particular, we have chosen to absorb 
the explicit powers of the lattice spacing a into the coefficients c, c. Consequently, as a function of a these coefficients 
scale according to 

cq,c ,c 3 ,c 3 ,c 5 ,c 5 = 0(a), c 2 ,c 2 ,c 4 = 0(a 2 ). (3) 
The quark mass matrix is given by 



m 



2m + m s 2(m-m s ) 



AI q = m \=M Q 1 + M 8 T H , M = — — , M s = -i 7 =^, (4) 



m s 



V3 



where isospin symmetry (m u = md = m) is assumed. Note that 0(a) contribution is already absorbed in the definition 
of M q , so there is no 0(a) term in the chiral Lagrangian. In the a — > limit, the pscudo scalar meson masses are 
related to m and m s according to 

{m\ = 2Bm, a = 1,2, 3, 

m% = B(m + m s ), a = 4, 5, 6, 7, ^ 
m 2 = (m + 2m s ), a = 8, 
3 

which, of course, agree with the continuum ChPT result. We also define the average mass 

1 x - „ 2m + m, 

— E< = 25 ^. (6) 



which will be a useful short hand notation in the following. 

Note that, except for a ir a independent constant, the term proportional to £5 is identical to the term with C3. 
Therefore we can set £5 = without loss of generality. 
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C. Shifted quark mass at leading order 



By expanding the LO Lagrangian to second order in 7r a , we obtain 

l {2) - 5 £KW (7) 

a 

The pseudo scalar meson masses at LO are therefore given by 

m 2 a = m 2 a (l-Nc 3 -£ 3 )-ml v Nc 3 -Nc 2 -£ 2 . (8) 

In the following we keep the number of flavors N arbitrary, but put N — 3 in the final expressions. 
We now define shifted quark masses, which satisfy 



2Bm, a = 1,2, 3, 



m 2 K = B(m + fh s ), a = 4, 5, 6, 7, 
2B 

m 2 = —^-{rh + 2m s ), a = 8. 



(9) 
(10) 



Explicitly they are given by 



ft at ~\ 2m + m s Nc 2 + c 2 , . 

m = m(l-Nc 3 -c 3 ) Nc 3 — — , (11) 

n at ~\ 2m + m s Nc 2 + c 2 

m s = m s (l -Nc 3 -c 3 ) A^c 3 — — . (12) 

o Id 

From these expressions the LO critical mass for N = 3 is defined by the condition 

m lv = m L( l ~ 6c 3 - c 3 ) - 3c 2 - c 2 = 0, (13) 

leading to 

2Bm cri ticai = 3c ^ + C2 _ = -Sml v . (14) 
1 - 6c 3 - c 3 

This definition for the critical quark mass assumes that all three quark masses are extrapolated to the massless point, 
including the strange quark mass. In numerical spectroscopy calculations, however, a different definition is sometimes 
employed where the strange quark mass is kept fixed at (approximately) its physical value. For this procedure the 
condition for the critical quark mass reads 

ml = ml (I - 3c 3 - c 3 ) - m^ v 3c 3 - 3c 2 - c 2 =0, (15) 

which results in 

on f \ 2Bm s c 3 + 3c 2 + d 2 

2Bra critica i(?7i s ) = . (16) 

1 - 5c 3 - c 3 



The difference between these two values is therefore 

/ \ c 3 
^critical ^critical \J^s J 



3c 2 + c 2 
1 - 6c 3 - c 3 



2Bm 6 



(17) 



1 - 5c 3 - c 3 

Indeed, numerical 2+1 flavor simulations |l5j | suggest a non-vanishing value for this difference with 

^critical - TOcritical("l s ) > 0. (18) 
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D. NLO Lagrangian 



The NLO Lagrangian provides the necessary counter terms in order to remove the UV divergences in the 1-loop 
integrals. The contribution of 0(M 2 ) is given by 



l nloSp = i4(S w )(M,S + E f M,) + L B (S W (M,E + ^M g )) + L 6 (M q Z + ^M q ) 2 
+ L 7 {M q Ys - Y^M q ) 2 + L 8 (M 9 EM g S + ^M q ^M q ), 



(19) 



where we introduced E^ = d^ES^Et and M q = BM q , where M q is the shifted quark mass matrix constructed from 
and rh s (cf. previous section) . The NLO constants Li are related with standard Gasser-Leutwyler coefficients Lf L 



atGL 



as L 4j5 = 2L^5 and ^6,7, 

The complete Lagrangian at 0(a 2 M) and 0(aM 2 ) is cumbersome and has not been computed so far. Here we 
only list the terms that contribute to the meson propagators, which we need for the calculation of the pseudo scalar 
masses. These terms are straightforwardly found by a spurion analysis with spurion fields proportional to the lattice 
spacing HHHH- Our result for the 0(a 2 M) La grangian reads 



L 



NLO, a 2 M 



Wi 

W 2 ' 



+ ^<E w (E + Et))(E + Et) 



2iV v 
f<E w (E 2 



(Et) 2 ; 



(Et) 2 )) 



+ W 5 {M q T, + T^M q ) + W 6 (Z + E t ) 2 (M g E + E+M g ) + VK 7 (E 2 + (E t ) 2 )(A/ g E + E+M g ) 
+ W 8 (E + Et)(A/ 9 E 2 + (Et) 2 A/ 9 ) + W 9 (M ? E 3 + (^) 3 M q ) + W W {Z + &){M q ) 
+ W n ((d^) 2 + (d^) 2 ) + W 12 ((d^) 2 (^) 2 + S 2 (5 M Et) 2 ). 



(20) 



while the terms of 0(aM 2 ) are given by 



7 NLO, aM 2 



= (E, 



2V (M q ) + j^(E + E+)(M 9 E + E+M g ) + y 2 (A/ g E 2 + (Et) 2 M g ) 



+ ^(E W (S + Et))(M ? E + EtM,) + ^(E^A^E + EtM ? ))(E + Et) 
+ F 5 (E, 1AI (M 9 E 2 + (Et) 2 Af 9 )) + 2V 6 {E^M q ) 

+ V 16 ((d^) 2 M q + M q {d^f) + 1/i 7 ((^E) 2 E t A/ 9 E t + EA/ g E(9 A1 E t ) 2 ) 



Et) 



y 7 (M|) + y 8 (A/ g E + EtM^) 2 + V 9 (M q T,M q E + YtM q YtM q ) 



+ (MqE + StAf g ) [y 10 (M g ) + Fn(A/ 9 E 2 + (Et) 2 A/ g )J + y 12 (A/ 2 E + EtM 2 ) 

+ V 13 (M q EM q E 2 + (Et) 2 M g Etiig + Vu(Z + E+>(M,E - EtAf g ) 2 
+ Vis^E 2 - (Et) 2 A/ 9 )(A/,E - E+M,). 



(21) 



As in the leading order Lagrangian we have absorbed the explicit powers of the lattice spacing in the low-energy 
constants, and their scaling behaviour is therefore 



Wi 
Vi 



0(a 2 ), 
0(a), 



= 0...12. 
0...17. 



(22) 
(23) 



III. RESULTS 



A. Quark mass dependence of meson masses 



The calculation of the pseudo scalar masses from the chiral Lagrangian in the previous section is straightforward. 
We collect some details and intermediate results of our calculation in the appendix. Here we simply quote the final 
result for the quark mass and lattice spacing dependence of the pseudo scalar meson masses. The total contribution 
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from LO tree, LO 1-loop plus NLO tree for and mx and are given as follows: 

m* )total = x + 2y+± [LUAZx + B^y + 5C} + V K {A\x + B\y + AC} + L r n {A*x + B«y + C} 
- {(Co + D )x + 2C y + (C w + D av )x 2 + D^y 2 + 2Dxy}] , 

X ~y+J2 [ L k{A%x + B*y + 6C} + Ll{A*x + B*y + 3C} + L r n {A*x + B«y + C} 
{(Co + Do)x - Coy + (C a v + D &v )x 2 + D*y 2 - Dxy}] , 

x - 2y + j 2 [L'{A*x + B*y + 3C} + L r v {A*x + B%y + 3C} + U K {A\x + B\y + AC} 



'if,total 



L rj, total 



- {(Co + Do)x - 2Coy + (C av + D av )x 2 + D^y 2 ~ 2Dxy}\ 



The parameters 



m 2 



— (2m + m s ), y 



IB 

—2BM 8 = —(m-m s ), 



represent the quark mass dependence. The chiral log is denoted by 



16tt 2 



V°g(™a), 



whose coefficients contain 



c = 








h 2C « 




a; = i(2c; - \x), 




K = 


\{2c : - 


In 


A\ = 


- 2X), 




K = 




5F), 


Bl = 




Y), B 


j = i(i2cB-n 




A K - 
B§ = 


>«- 


3X), 
+ 3Y), 


A« = 

7T 




-§*), 




6c 5 ), 


K = 


\(^- 


2 h 


Al = 


\{^- 


- 3 -X), 
2 h 


A\ = \{2Cl-2X), 




B« = 






24c 5 ), 


Bl = 


i(-3F- 


f36c 5 ), Bl = ±(-3Cl 


f 5Y 



These terms are parametrized by 



X 

ci 



L n 



A 



1 - 6c 3 - 45 3 - 36cjB , Y = (1 - 3c 3 - 45 3 )£, Z = [9c 2 + 45 2 + <5m^ v (l - 18c 3 - 4c 3 )] , 



2+18c + 9S + 15c 3 B, C£ 



1 + 24c + 65 + 12c 3 B, C™ = 6c + 5 + 3c 3 B, 



K 



= - + 24c + 9co + 18c 3 B, C 
= 6c + 35 + 9c 3 B, C'l = 18c 



9. 



7 + 18c + + 9c 3 B, C* = j + 6c + ^cq + 3c 3 B 



with 



A = 



1 



C3 



4 
3c 



S = 



9c 3 B, C\ = 3 + 24c + 10c + 12c 3 B, 



1 



1 - Nc 3 ~ h 



1 - 2AT C3 

The polynomial (non-logarithmic) terms contain 

16L 5 + 8(L 8 + 4)+48Va, D«=4L 5 
16L 5 + 24Z 8 + 24L^ + 96Z7 + 48y A , D 



Sm axr = -■ 



Ncn 



C2 



D 7T 

yy 
yy 



C av = A(NL 6 



C 



NL 4 



L 5 ), Aw = 2(i 8 +L 8 ; 



20L 8 

- C av 



1 - 27Vc 3 - c 3 

- 16L 8 + 48Va, 
AL 5 + A(L S + L' 8 ), 



o — 



D 



A(W + W x + W 2 + W 3 + W 4 ) + 2W 5 + 8N 2 W 6 + ANW 7 + 16NW S + 18W 9 

~8(W 11 + W 12 )^a 2 W c , 

N[16(NW 6 + W 7 ) + AW S + 2W W ] = a 2 W D , 



where 

Li = Li + V + Vi + V 2 + V 3 = L A + aL\, L 5 = L 5 + V4 + V 5 + V 6 - Vie - V17 = L 5 + aL\, 
L e = L 6 + 2NV S + ±V W + |fu = L 6 + aL\, L 7 = L 7 + 2NV 1A + 2V 15 = L 7 + aL\, 

L 8 = L 8 + 2iVy 9 + iy 12 + |fi 3 = L 8 + aLl a , L' s = L 8 + 2NV 9 + 2V 13 = L 8 + aLf , 

Kv = N(V r + 2Vg + W 8 ), V A = ^(V 7 + 2V 9 ). 

Even though the NLO parameters have been used to remove the divergent terms from the loop integrals we use the 
same notation for these coefficients. We finally note that there exists the following constraint among some of the 
coefficients: 

J2 A b = ;r[6 + 96co + 32co + 60c 3 5-5X]. (33) 
fa 3 

Consequently, in the limit y — > we obtain identical results for m„ or m„, as it should be for three degenerate 
quark masses. 

Obviously the final results for the pseudo scalar masses are fairly lengthy. From a practical point of view the number 
of independent unknown parameters in these expressions is crucial for their usefulness. Unknown parameters arc the 
critical quark mass m C riticai, the constant 2B and the decay constant /. The coefficients of the chiral log terms, given 
in ea. (|29|l - H32|l . are defined through five independent 0(a) parameters, cq, 60, C3, c 3 and C5, and the coefficient C, 
which is an independent 0(a 2 ) parameter. In the analytic NLO correction we find the independent combinations 
C av + D av , Dy y , Dy y , D ^ and D, which start at 0(1), and the two coefficients Co, -Do are of 0(a 2 ). Therefore, the 
total number of independent parameters is thirteen besides m cr iticai) 25 and /. 



B. 0(Ma) term at NLO 



In eqs. (|24fl . (|25|l and i|26|) . 0(Ma) terms are considered at LO. In this subsection we present the results for the case 
that these terms are treated as NLO corrections. First of all, the expressions for the shifted quark masses simplify to 



Nc 2 + 5 2 
m — — = m - mcriticai) 



2B 
Nc 2 + c 2 
2B 



ill's '^critical ' 



The quark mass dependence of the pseudo scalar meson masses becomes 



where 



4,totai = (x + 2y)[l - Nc 3 - c 3 ] - xNc 3 + 



1 



Ll{-x + y + 5C} + U K AC + LU-x - -y + C} 



{(Co + D Q )x + 2C y + (C av + D av )x 2 + D; y y 2 + 2Dxy}] 



'if.total 



= ( x -y)[l- Nc 3 - c 3 ] - xNc 3 + — 

f 

{(Co + D )x - C y + (C av + D av )x 2 + D*y 2 - Dxy}] , 



U K 6C + Ll{-\x + SC} + n v {\(x -y) + C] 



<totai = (x-2y)[l-Nc 3 -c 3 \-xNc 3 + — 



L r v {~\x + \y + 3C} + Ll{-\x -y + 3C} 



+ L r K {-(x -y) + AC} - {(Co + Do)x - 2C y + (C\ v + D av )x 2 + D^y 2 - 2Dxy} 



C 



Z _ (9 - N)c 2 + 3c 2 
6 ~ 6 



yy 



1 6L 5 + I6L 8 = AD 



K 

utr 



D« y = 16L 5 + 48L 8 + 96L 7 , D = C av + AL 5 + 8L 8 , 

C av = A(NL 6 + NL 4 + L 5 ), D av = AL s , 

C = A(W + W 1 +W 2 + W 3 + W A ) + 2W b + 8N 2 W 6 + ANW 7 + 16NW S + 18W 9 

- 8(W n +W12) = a 2 W c , 

Da = N[16(NW 6 + W 7 ) + AW S + 2W W ] = a 2 W D . 



(34) 
(35) 



(36) 



(37) 



(38) 
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The number of independent parameters is reduced compared to the result given in the previous section. Besides 
^critical) 2-B and / there are C3, £3, C, L4 + Lq, L$ + Lg (note that C av + D av , D yy , D yy and D can be expressed by 
L 4 + L 6 and L 5 + Lg.) and D^ y . The total number of independent parameters besides m cr iticai, 25 and / is reduced 
from 13 to 6. 

However, for improved theories there is some lattice spacing dependence implicit in the definition of the renormalized 
quark mass. This results in 3 parameters (b B + bm , bm , bm ), as we will show in the next subsection, where we 
discuss 0(a) improved theories. 



C. Formula in O(a) improved theories 



We finally consider the case that a non-perturbatively 0(a) improved quark action (i.e. the clover quark action ) 
is used in the lattice simulations [Til UtI ITa. Il9l l2pj . In this case there are no on-shell 0(a) terms in the Symanzik's 
effective theory provided that the relevant improvement coefficients are tuned non-perturbatively to an appropriate 
value. In particular, 0(a) improvement requires that some a dependence is absorbed in the definition of renormalized 
masses and the gauge coupling: 



m 

m., 



m + o m Wa + b$(2m + m s )ma + b m ^(2m 2 + m 2 s )a, 
m s + b$m 2 s a + b$ (2m + m s )m s a + b$ (2m 2 + m 2 )a, 
(2m + m s )a N 



9l 



(39) 
(40) 

(41) 



where b g and b 



3(fr 



(2) 



6m') are improvement coefficients defined in Ref. [19|,|20|. Therefore, as long as 
on-shell quantities are considered, there are no terms of 0(a), O(Afa), 0(Af 2 a) etc. in the WChPT Lagrangian, if we 
replace 



fh — > 77i = m + b$m 2 a + i)j^'(2m + m s )ma + &^'(2m 2 + m 2 )a, 

fh s = fhg + b^7n 2 a + b$ (2m + m s )m s a + b$ (2m 2 + m 2 )a 



m s 
D 



B = B[l + b B (2m + rh s )a] 



(42) 
(43) 
(44) 



Here the last modification comes from the mass dependence of <? 2 in the Symanzik's effective theory. 

We emphasize that there are no terms linear in a in the chiral Lagrangian and in the results for the pseudo 
scalar masses as long as one parameterizes it in terms of m, which absorbs some 0(a) dependence through proper 
renormalization. Using 777 instead, which is simpler in practice since this mass is directly proportional to the difference 
between the bare and the critical quark mass, there is some 0(a) dependence left explicit. 

Having made these remarks we can write down the WChPT expressions for non-perturbatively 0(a) improved 
theories: 



m 7T, total 



TO _R",total 



2y + j- 2 



+ y + 5C} + U K AC + L^{ l -x -\y + C] 



{(C + D )x + 2C y 
{(C + D )x - C y + (C* av + D av )x 2 + 



+ (C av + D av )x 2 + D; y y 2 + 2Dxy}] , 
1 
3 

Dxy}] 



Ll{--x + 3C} + L r v {^x-y) + C} 



<y 2 



m ri, total — x 



l 

P 



L r v {- 



1 



1 



-,y - 



3C} + Ll{--x 



{(Co + Do)x ~ 2C y + (C a v + D„)x* + D^ y y 2 



2/ + 3C}4 
- 2Dxy}] , 



y)+4C} 



(45) 



(46) 



(47) 



where 



2B /n 
(2m 

2B 



m s ) = x[l + (b B + b$)(2m + m s )a] + (b% + 36 m ))(2m 2 + m 2 )a, 



(m — m s ) = y[l + (b B + b l ^ ) )(2rh + rh s )a + b$(m + m s )a] 



(48) 
(49) 
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and 



c 


Z (9-N)c 2 + 3c 2 _ 
= 6 = 6 > ^ = 16L * + 


16L 8 = 4L>f y , 


D v 

yy 


= 16L 5 + 48L 8 + 96L 7 , D = C av + 4L 5 


+ 8L 8 , 


C av 


= 4(NL 6 + NL A + L 5 ), D &V =4L 8 , 




C 


= 4(1U + IUi + TU 2 + 1U 3 + Wi) + 2W 5 4 
- 8{W n + Wi 2 ) = a 2 IU c , 


- 8N 2 W 6 + 4AIU 7 + 16AIU 8 + 


D 


= N[16(NW 6 + 1U 7 ) + 4W 8 + 2Wio] = a 2 





Independent parameters (besides ^critical, 2S, and /) are 65 + 6; n . 6m , bm , C. L4 + Le 5 -^5 + -^8 an d -D^- The 
number of independent parameters besides m cr iy ca i, 2B and / is therefore 7, reduced from previously found 13. 

IV. CONCLUDING REMARKS 

In this paper we computed the pseudo scalar masses in 2+1 flavor WChPT. We presented results for three different 
order countings, appropriate for various sizes of the O(aAf) term in the chiral Lagrangian. Depending on the lattice 
action used in the numerical simulation (unimproved, pcrturbatively improved, non-perturbativcly improved) one has 
to choose one result for the chiral extrapolation. Since we have no prior knowledge about the size of the 0(aM) 
contribution we suggest to perform fits to the data with all three forms and let the data decide which form is most 
appropriate. 

The number of unknown fit parameters is significantly larger than in 2 flavor WChPT. Using our results requires 
sufficiently enough data points in order to perform the chiral fits. The CP-PACS/JLQCD collaboration is currently 
performing 2+1 flavor simulations at three lattice spacings with five values for the light up and down quark mass and 
two different strange quark masses. At least for these simulations the number of data points exceeds the number of 
unknown fit parameters. Performing the chiral extrapolation of the CP-PACS/JLQCD data using our results is work 
in progress. 
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APPENDIX A: USEFUL FORMULAE 
1. Expansions in terms of n 

In this subsection we collect some useful formulae necessary for the expansion of the chiral Lagrangian in terms of 
the 7r fields. 
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a. LO terms 

At LO we have to expand terms at 0(7r 4 ). 

4 8 
(d^Ed^E 1 ) = —{d^d^) + — jid^nln^d^Trln), 

(ME + £tM) = (2M) - j 2 {M^) + ^(Af^ 4 ), (S + Et) = (2) - + A (tt 4 ), 

<S + Et) 2 = (2) 2 - A<2)<^) + ^(2)(. 4 ) + ±J(. 2 ) 2 , ((E + Et) 2 ) = (4) - ^(n 2 ) + ^(tt 4 ), 



(E + Et - 2)(a Al E3 Al E t ) = -^{^(d^dpiT), ((E + Et - 2)^E9 M Et) = -H^^^), 



(E + Et) (ME + EtM) = (2)(2M) - -|(M)(tt 2 ) - - (1)(Mtt 2 ) 



/ 2N " ' / 2 
16 

T 4 ^" 1 ' 3/ 4X "" ,v ' ' ' 3/ 4 



+ 7i (7r 2 )(M^ 2 ) + -j(M)(^ 4 ) + — (1)(A^ 4 



((S + Et)(MS + EtM)) = (4M) - ^(Mtt 2 ) + ^(Mtt 4 ), 
<E - Et) 2 = 0, ((E - Et) 2 ) = -i|(. 2 ) + J* 

(E - Et) (ME - EtM) = ^(n 3 )(Mn), ((E - Et)(ME - E*M)) = -^(Mtt 2 ) + ^(M^). 



b. NLO terms 



We have to expand the NLO terms to 0(7r 2 ). 



(E w )(ME + EtM) = j^{d ll nd fl n){2M), (E W (ME + E f M)) = — {d^d^M), 



8 16 

(ME + E f M) 2 = (2M) 2 - — (2M}(M7r 2 ), (ME - E f A/) 2 = - — (Mtt) 2 , 

8 16 

(MEME + EtMEtM) = - — (MttMtt + MV), (E 2 + (£ f ) 2 ) = (2) - — (tt 2 ), 

/ / 

(ME 2 + (Et) 2 M) = (2M)-^(Mtt 2 ), (ME 3 + (E f ) 3 M) = (2M) - ^(Mtt 2 ), 

J J 

(M 2 E + EtM 2 ) = (2M 2 ) - -^(MV), (MEME 2 + (S f ) 2 MEtA/) = (2M 2 ) - 4r(4M7rM7r + 5M 2 7r 2 ), 

J / 

(ME -E f M) = j(Mtt), (ME 2 - (Et) 2 Af) = j(Mtt), 

((d^r + (a^Et) 2 ) = {{d^f{^f + T 2 {d^f) = -^{d^d^). 

2. Formula for traces 

After expanding the Lagrangian in terms of the 7r fields, we have to take the trace in the flavor space. 
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a. LO terms 



(1) 



2B(M q ir 2 



2B(M q TT 4 ) = 



abed 



^ F abcd 7T a 7T b 7T c 7T d . 

a,b,c,d 

2BM | 2 





^ — ^ 2 

a 






[ i 


= 2Bm, 




= 1,2,3, 


1 m 2 K 
{ m v 


= B(m + m s ), 
2B , 

= — [m + 2m s ), 


« 

a 


= 4,5,6,7, 




1 \ ^ 2 2 


1 

f 8 


d abe d cde ir a iri,ir c ird 



.b 





a~e a,b,c 


afcc 7r a ^7r c , 2B(M g7 r) = ^2BM 8 n s 


NLO terms 






(2BM q ) 


= 7Vm 2 v , (^) = -^7^, 
a 




{2BM q ) 2 


= iV 2 (m 2 v ) 2 , ((9^) 2 2BM 9 ) = m 2 y i E^^) 2 

a 


+ ^ m 2 \Y,d aa \d^ a f, 

a 


((2BM q ) 2 ) 


= 7V(m 2 v ) 2 + i(Am 2 ) 2 , (2BAV 2 ) = ^E^ 

a 


2 , (2BAI q ir) 2 = i(Am 2 ) 2 7r 8 2 , 


2((2BM q ir) 2 } 


= (mL) 2 E^ + -LAm 2 ^^ 2 + i--p! 

a a 


|- 8 2 +E{( daa8 ) 2 -E(/ Qb8 ) 2 }- 

a b 


{(2BM q ) 2 ir 2 ) 


- ^(-L) 2 E^ + f» A ™ 2 E r8 ^ + (a- 2 ) 2 E - 

a a a 



3. Group factors 

We have to calculate some Lie group factors. 

(ja>b — ^ jabc jabc ^ S^^C^^ C^ a 

c 

C™ = C KlKl - C K2K2 = 1, C nK = C KlK2 = 1/4, C nv = 0, C Kr] = 3/4, 
where 7r represents a = 1, 2, 3, K represents a = 4, 5 (i^i) and a = 6, 7 (-K2) 77 represents a = 8. 

jj&b — ^ ^aac ^bbc ^jaa3 ^jbb3 _|_ ^aaS ^bb8 £)^ a 

c 

L)^ = D KlKl = D K2K2 = D m = 1/3, D nK = D KlKo - = -1/6, = -1/3, = 1/6. 

L> Qb = J2d abc d abc = Z) ba , 

C 

D aa = 1/3, 5^ = £> A ' lA ' 2 = 1/4, 0^ = 1/3, D A '' 7 = 1/12, others = 0. 
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E ab = V d ccS d aac d bbc = -^=\d aa3 d bb3 ~ d aa8 d bbs ] = E ba , 

J^ktt _ j^l<il<2 _ ^7"?f; _ __}_ JJJ"K j^KiKi j^K 2 K 2 1 A 7 ™' = - A^"'' - 

3^3' " ~ 6-n/3' ~ 3^' 6^' 



E ab = ^d ccB d abc d abc = E ba , E aa = -== <| 1/6 for A', 



-, [ —1/3 for 7r, 
" /6 for A'. 
1/3 for r\, 



1 



A 7 '* = i=, A 7 '" 1 ' = -i=, A^ 2 = A K " = L,, others 

8^3 3\/3 4V3 24V3 



= 0. 



With these definition we have 



4i? 



aabb 



2BM {— + -D ab 
^ N 2 



2£M <j + Il) Qf > 



The following formulae for TV = 3 arc useful. 



25M 8 <! i_<5 o6 d aa8 + \e° 



N 



2N 

where a = 1,4, 



-, -, / 5/6 1/2 1/2 1/2 1/2 1/2 1/2 1/2 \ 

— + -A> ab + A> ab = 1/2 1/2 1/2 5/6 1/2 1/2 1/2 1/2 , 

N 2 \ 1/2 1/2 1/2 1/2 1/2 1/2 1/2 5/6 / 

/ 1/3 1 1 1/4 1/4 1/4 1/4 

(d aaS + d m ) + -E ab + -E ab = 



_ I 1/4 1/4 1/4 -1/6 -1/2 -1/2 -1/2 



5/4 I , 

\ 1 1 1 -5/4 -5/4 -5/4 -5/4 -7/3 



4. 1-loop contractions 

The following contraction formula is useful for the calculation of the meson propagators at 1-loop. 
(d,n[n,d,7T}n) - -J]TC ab [(^7r„) 2 7 (mg) + ^(m 2 )] , C ab = £/ abc r bc , 



/ (m 2 ) 



16tt 2 



[A + l-logm 2 ], A = 7 + log(47r), A(m 2 ) = -m 2 I (m 2 ), 



(2BMn i ) -> ^7r22/ (m 2 )[A oabb + 2A 



ababi 



a . & 



™2\ 



^(1 + 2<5 afc ) + i ^(d QQe d bbe + 2d abe d abe ) 



\Y J *lh{rhl){\ + 25 ab ), 



a . b 



(n 2 )(2BMn 2 ) -> I £ 7r 2 / (m 2 )(™ 2 + m 2 )(l + 2J a6 ), 

a, 6 

- jE[(^a) 2 /o(m 2 ) + 7 r 2 A(m 2 )], 



a . b 



4iV 8 



(ir 3 )(2BMir) -2AM 8 d aa8 [/ (m 2 )^ 8 2 + I (m> 2 ] 
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e - J cfixLiofro) I T , 7r ^ p -/d 4 xL( 2 '( 7rQ ) p -/d 4 2; (L< 4 )( 7rQ +7r G )-L( 4 '( 7 r G )) 



S eff (7r G ) = / d 4 xL LO (7r G )-log( e -^ 4 ^ i(4, ^ + ^)- i(4> ^»), (B4) 



APPENDIX B: CALCULATION OF THE PSEUDO SCALAR MESON PROPAGATOR AT 1-LOOP 

In this appendix we give some details of the calculation of the pseudo scalar meson propagator at 1-loop in the 
WChPT. 

1. Effective action for BG fields 

In order to calculate meson masses at 1-loop, we use the background (BG) field method. Wc first split the tt field 

as 

7T = TTQ + 7T G , (Bl) 

where ttq represents quantum field while ttq is a BG field which satisfies the equation of motion. Inserting this into 
the chiral effective Lagrangian. wc have 

£loM = iLo(7r G ) + L (2) (^Q) + i (4) (^Q+^G)-i (4) ^G). (B2) 
Integrating out ttq wc obtain the following formula 



This leads to 



where 

</(tq,tg)> = Jv7r Q e-f dixLl2) ^f(Tr Q ,7r G ). (B5) 
Expanding S e & in terms of 7Tg, wc obtain 

oo 

S eS (tt g ) = constant + S$} (tt g ) + ^ sjg (tt g ) , (B6) 

n— 3 

where S^g (ttg) contains the n-th power of the field 7r G . In the calculation of the pseudo scalar masses wc arc interested 
in the n = 2 case. We write 

S&VeO = / d 4 xL ho (n G ) + s[% oop (7r G ) + ■■■, (B7) 

where • ■ • represent the higher loop contributions. We call Sj__ lo the 1-loop contribution to the meson propagator 
and write 

s[% oop («g) = J d*xL<?l loop (v G ). (B8) 
2. Expansion of the LO Lagrangian 

We now expand the LO Lagrangian Q in terms of the pseudo scalar field 7r a . Using the expansion and trace 
formulae given in appendix ^ we obtain 

= {{d^f) +2B{Aiy) -|(2)<^) -c- 2 (7r 2 ) 

- c 3 (2BM q )(TT 2 ) - c 3 {l)2B(M q n 2 ) - c 3 2B{M q -K 2 ) 

= ^E^) 2 +m 2 7r 2 ], (B9) 

a 
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at second order in 7r a , where the pseudo scalar meson masses at LO are given by 

ml = m 2 a {l-Nc z -c 3 )-m 2 m Ncz-Nc 2 -c 2 , (BIO) 



N 2 



For N = 3 flavors we have 



= 2Bm, a = 1,2, 3, 

4 = B(m + m s ), a = 4, 5, 6, 7, 2 _ 2m + m s 



< = < »*jr-^"*-r»*.;, i*-*,o,o,i, m g v ^ 2 g ^ — , (B12) 
m 2 = -^-( m + 2m s ), a = 8. 



Ea. (|B9|) gives the pseudo scalar meson propagator at LO. 
The 4-th order terms in the LO Lagrangian become 



+ ^ + + NC3mlv ^ + + ^V)(2^ 2 ) 

- ^VX(^) 2 ) - ^V(M 2 ) + ^^)(2BMir). (B13) 
There terms give the 4-point interaction vertices of the pseudo scalar mesons. 

3. 1-loop contribution to the propagator 

Using the formulae in appendix [S] it is now easy to calculate L^]_y . Including the tree level contribution we 
obtain 

S { c 2 s \n) = J d 4 4]Tz a [(^ a ) 2 + miX] . (B14) 

a 

For the 7T (a = 1, 2, 3), K (a = 4, 5, 6, 7) and r\ (a — 8) we find the wave function renormalization as 



Z v = 1-— j [L ff {2 + 9(2c + oo)} + ijr{l + 6(4co + ao)} + i,(6co + oo)] 1 

Z K = 1 - i [ix{3/4 + 9(2c + c Q /2)} + L K {3/2 + 3(8c + 35 )} + L„(3/4 + 6c + 5c /2)] , 
"J 

Z„ = 1-— J [i w 3(6co + c )+iK{3 + 2(12co + 5co)}+i, ) 3(2co + co)], (B15) 



where 

Ltt = Io(ml), L K = h(m 2 K ), L v = h{m 2 .). (B16) 
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Similarly we have 



ml 



m K,R 



' K 



L r),R 



ml 



1 

w 



o 5 o Am 
Cllml - -Xml - 5Y^= 



+ L K {Cl{m z K + mi)-2Xm 



"^av 



T,2 



Y 



.Am 

i7f 



2Z 



Am 



7,2 



+ L, ^ C^m* + ml) - :im^ — (Y — 12c 5 )^j= + \z 



1 

3p 



4-y/3 



3Z 



^2 
av 



-(y — c 5 ) — =■ + -z 

4 V 5 4V3 2 



1 



< Cp.m% - Jlmi + 5(Y - ^c 5 ) 



f 3 Am 2 3 1 

L w |^(m 2 + ml) -Xml 3(K - 12c 5 )^ + ^ j 



4^/3 2 
4V3 ' 2 



24 

L x <{ ^(m 2 - + m 2 ) - 2Xm 2 v + 5(Y - — c 5 



2Z 



(B17) 



(B18) 



(B19) 



5 4^3 

for the pseudo scalar meson mass. The parameters in these expressions have already been given in subsection IIII Al 



4. NLO contribution to meson propagators 

Using the formulae for the expansion in powers of the pion field and the trace formulae in the appendix^ we have 



iNLO = -C°[(3^a) 2 +™l NL o^] 



(B20) 



where the wave function renormalization factor is given by 



^NLO 

'a 

'"A 



-p [™av z av + Am 2 Z^ + Zq] , 

4(7VZ 4 + i 5 ), 
2d aa8 L 5 , 

4(W + Wi + W 2 + W 3 + W A ) - 8(W n + W 12 ), 



and the mass term is defined as 
1 



m o,NLO 



2 [m 2 a {m 2 v Cav + Am 2 C a A + C } + m 2 w {m 2 v C av + Am 2 D a A + D } + {Am 2 ) 2 E a A ] , (B21) 
4{NL 6 + NL 4 + L 5 ), D av = 2(L s + L' s ) + V avi C a A = 2d aa8 L 5l 
2d aa8 (L 8 + Lg), E a A = e a L s + e' a L' s + V A + 5 a8 f.2L 7 + ±L' S , 



4{W + ^1+^2 + ^3 + Wi) + 2W 5 + 8N 2 W 6 + ANW 7 + 16NW 8 + 18W g 
-8{W ll + W 12 ) = a 2 W c , 



Co 

D = 2V[16(AW 6 + W 7 ) + AW 8 + 2W W ] = ^W, 



1G 



The constants here are given by 



and 



(B22) 
(B23) 

(B24) 

(B25) 

(B26) 

(B27) 



V m = N(V 7 + 2V 9 + iNVs), V A = hv 7 + 2Vg), (B28) 

e« = ^-^ daa8 > (e„e K ,e n ) = (l, A (B29) 

e' a = \[(d aaS ) 2 -£(/ a " 8 ) 2 ]> {e'«,e' K ,e' ri ) = {\-\, 1 -). (B30) 

b 



U 




V V + Vt 


+ v 2 + 


^3 = ^4 + 0-^4, 


f 

Li 5 




h V4 T v5 


J- T/ 
+ V6 ~ 


Vl6 - Vi 7 = 


r i r 1 

As + aL 5 


U 


= L 6 - 


h 2NV 8 + 


1 


- ^ii = L B 


+ aLl, 


L 7 


= L 7 - 


h 27Wi 4 - 


- 2V 16 = 


= L 7 + aL 7 , 




L s 


= L 8 - 


h 2W 9 + 




5 

' ~Vl3 = ^8 


+ aLl a , 


L' 8 


= Ls- 


1- 2NVg + 


2^13 = 


L 8 + aLg 6 , 





5. Cancellation of UV divergence 

In order to perform the renormalization, we consider divergent parts of meson masses in eq. I|B17|) - l|B19|) . which 
are given by 

Kkv. = - [ C > + C > + C ^ 2 + C yyV 2 + C *v x v] > ( B31 ) 

where x = m^ v and y = j^Afn 2 , in terms of which, = x + 2y, fh 2 K = x — y, and = a; — 2y. Constants are 
given by 

ci = cf = c;> = 5z, c™ = 2z, = -z, q; = -2Z, 



= 2^C 6 Q a -5X = 2(3 + 48c o + 165o + 30c3B)-5X, 



C£, = 8CZ - CI - 7Y - 2Ac 5 B = 15 + 120c + 66£ + B[-7 + 129c 3 + 285 3 - 24c 5 ], 

= 2Cf + 2Cf + 6C* - 7y + 12c 5 S = 9 + 120c + 425 + B[-7 + 93c 3 + 28c 3 + 12c 5 ], 
= &C% + 3C£ - 21Y + U4c 5 B = 9 + 120c + 545 + B[-21 + 171c 3 + 845 3 + 144 Cj 



5 J 



C£„ = 8C; - C£ - 4C,^ - 2X - 7Y + 12c 5 B = 15 + 96c + 62S + 96c 3j B - 2X - 7Y + 12c 5 B, 

7 - 1 

-y- 6c 5 5 = --< 
2 2 

= -8C* + 4C% - bCl + 2X + 7Y- !2c,B = -<%,. 



Cg = -4Cf + 5Cf - 7C* + X + -Y - 6c 5 B = --<%,, 



On the other hand, the NLO contributions lead to 



klNLO = ~ [D a x x + D a y y + D a xx x 2 + D a yy y 2 + D a xy xy] 



where 



D a x = C + D , Dy = 2C , D* = -C , D% = -2C , D a xx = C av + D 



av • 



Dy y = 16i 5 + 8(i 8 + 4)+48y A , £>^=4L 5 + 20L 8 -16L 8 + 48F A , 
Dly = WL 5 + 24L 8 + 24L' 8 + 96L 7 + 48V A , 
Dl v - 2C av + 8L 5 + 8(Z 8 + L' 8 ) = -2D% = -D^. 
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In order to cancel the UV divergences, the divergent part in the NLO terms must be chosen according to 



[ C o] div . 

[Cav + £>av] div . 

16L 5 + 8(L 8 + L' s ) + 48 Va 

- div 

4L 5 + 20L 8 - 16L' S + 48Va 

- div 

16L 5 + 96L 7 + 24(L 8 + L' a ) + 48 Va 

- div 



A 

"48tt 2 
A 



Z, [D ] c 



A 
'48tt 2 



r4Z, 



48tt 2 

A 
48^ 

A 



[2(3 + 48c + 165 + 30c 3 5) - 5X], 
2 [15 + 120c + 665 + B(-7 + 129c 3 + 28£ 3 - 24c 5 )], 



48tt 2 

A 

"48^ 
A 



[9 + 120c + 425 + B(-7 + 93c 3 + 28c 3 + 12c 5 )], 
2 [9 + 120c + 545 + B(-2l + 171c 3 + 84£ 3 + 144c 5 )], 
[15 + 96c + 625 + 96c 3 S - 2X - 7Y + 12c 5 B}. 



div. 487T 2 

Notice that we can remove all divergences [m^]div. consistently by these parameters, as it should be. 
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